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ABSTRACT 



Rate coefficients for inelastic Mg+H collisions are calculated for all transitions between the lowest seven levels and the ionic state 
(charge transfer), namely Mg(3s 2 >S, 3s3p 3 P, 3s3p 'P, 3s4s 3 S, 3s4s 'S, 3s3d 'D, 3s4p 3 P)+H(ls) and Mg + (3.s 2 S)+FT. The rate 
coefficients are based on cross-sections from full quantum scattering calculations, which are themselves based on detailed quantum 
chemical calculations for the MgH molecule. The data are needed for non-LTE applications in cool astrophysical environments, 
especially cool stellar atmospheres, and are presented for a temperature range of 500-8000 K. From consideration of the sensitivity 
of the cross-sections to various uncertainties in the calculations, most importantly input quantum chemical data and the numerical 
accuracy of the scattering calculations, a measure of the possible uncertainties in the rate coefficients is estimated. 
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1. Introduction 

The need for accurate data on inelastic atomic collisions in- 
volving hydrogen atoms for non-LTE line formation calcula- 
tions in cool star atmosphere s has been well documented (e.g . 



tions in cool star atmospneres nas been well documented (e.g. 
Steenbock & Holweeeii fl984t iLambertl [1991 iKiselmanl l200lt 



Asplundl I2005I) . Non-LTE models have often, for lack of a 
better alternati ve, been forced to use t he so - called "Drawin 
formula" (see ISteenbock & Holwegerl 1 1984 iDrawinl Il968t 
iDrawin & Emardl|l973l) . which has been shown to not contain 
the relevant physics nor provide order-of -magnitude estimates 
of the collision rates (Barkle m et al.ll201lb . 

Full quantum scattering calculations based on quantum 
chem ical data have been done for Li+ H dBelvaev & Barklerrj 
l2003h and Na+H dBelvaev et all l2010h . and rate coefficients 
provi ded for astrophysical applications dBarklem et al.l [2003, 
2010 ). The da t a have bee n applied to non-L TE calculations by 
Bark lem et al.l(l2003h and ILind et all (120091) in the case of Li, 
and ILind et all d2011l) in the case of Na. In both cases, it has 
been found that direct excitation and deexcitation processes, 
X(n/) + H ^ X(n'l') + H where X is the atom of interest, are 
of little consequence for the statistical equilibrium. However, 
charge transfer processes, X(n/) + H ^ X + + H~, particularly 
where nl corresponds to the first excited S -state, are important. 
The resulting effects on abundance corrections in solar-type stars 
were often small, of order a few 0.01 dex. In giants the influence 
was larger, of order 0.05 dex, and even larger still in metal- poor 
stars reaching of order 0.1-0.2 dex. Thus, these processes are 
important for accurate absolute abundances, and also for accu- 
rate relative abundances among dissimilar stars. 

Alkalis have naturally been the first cases studied since they 
have only a single valence electron and thus present the easiest 
cases for quantum chemical calculations. Also, Li and Na are of 



significant astrophysical interest. We now present data for Mg, 
the motivation for which is two-fold. First, it is of interest to 
study hydrogen collisions in more complex atoms and investi- 
gate the similarities and differences to the cases of alkalis: an 
atom with two valence electrons is the next logical step. Second, 
Mg is of significant astrophysical importance. Mg creates some 
of the strongest absorption lines in stellar spectra and is rela- 
tively easily detected even in low-quality spectra or in metal- 
poor stars. It is a tracer of the star formation history and chemical 
evolution of the galaxy since it is thought to be formed predom- 
inantly in supernovae resulting from massive stars. In fact, Fe is 
often used as reference element in galactic chemical evolution 
studies since strong lines of Fe are ubiquitous in the visual spec- 
tra of cool stars. However, a number of different astrophysical 
sites are expected to have contributed significantly to the pro- 
duction of Fe and thus its history is quite complicated. Given 
its relatively simple formation history, Mg is now a comm only 
used alternative reference element (e.g. ICavrel et al.l l2004). We 
note that, to our knowledge, there are no previous theoretical or 
experimental studies of Mg+H inelastic collisions. 



Quantum chemical calculations of the potentials and cou- 
plings have now b een done for MgH including excited states 
(iGuitou et al. 2010). These were then used in full quantum scat- 
tering calcu lations, initially testing on a small number of states 
(IGuitou et al. 2011), and later applying them to t he seven low- 
est sta tes of symmetry 2 2 + plus the ionic channel (Be lvaev et al.l 
|2012|) . This allows cross-sections to be calculated for direct ex- 
citation between the seven lowest states of Mg and for charge 
transfer processes involving these seven states. These calcula- 
tions should be sufficient for non-LTE line formation modelling 
of spectral lines corresponding to transitions between low-lying 
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levels in cool stellar atmospheres. Such models require rate co- 
efficients and these are presented here. 

2. Collision rate coefficients 

The rate coefficients, (crv), for excitation and deexcitation pro- 
cesses, 



Mg(3snl z " +l L) + H(ls)^Mg(3sn'l' z " +1 L') + H(ls), (1) 

and for the charge transfer processes, ion-pair production and 
mutual-neutralisation, involving the ionic state, 



Mg(3s nl zb+l L) + H(ls) ^± Mg + (3s Z S) + H , 



(2) 



are presented in Table Q] The coefficients have been obtained by 
folding the cross-sections with a Maxwellian velocity distribu- 
tion from thresholds to 17 eV in total energy (corresponding to 
at least 10 eV in collision energy), and are presented for temper- 
atures in the range 500-8000 K. 

The accuracy of the rate coefficients is determined by the 
ac curacy of the cross-se ctions near the threshold. As discussed 
in Barkler n et alj ([20 10), this accuracy depends on a number of 
aspects, most importantly the accuracy of the input quantum 
chemical data (potentials and couplings), the possible contri- 
bution of neglected molecular symmetries and states, and the 
numerical accuracy of quantum scattering calculations. As in 
Barkl em et al.l (1201 Oh . we have considered these points in detail 
for each transition and attempted to estimate"fiuctuation factors" 
for each transition, which should be indicators of the uncertainty 
in the rate coefficients. The probable range for a given rate co- 
efficient can be estimated to be from the value given in Table Q] 
multiplied by the minimum value of the fluctuation factor, up to 
the same value multiplied by the maximum value of the fluctua- 
tion factor. Table [2]presents minimum and maximum fluctuation 
factors for the rate coefficients in Table Q] 

The main considerations used to determine these fluctuation 
factors are as follows. Based on comparison of calculations us- 
ing different quantum chemical data sets from different calcu- 
lations done by us, we estimate the error from this source at 
of order 10%. Calculations were also done using two different 
programs for the quantum scattering calculations and the differ- 
ences did not exceed 10%. Together, this results in minimum and 
maximum fluctuation factors of 0.8 and 1.2, respectively, for the 
majority of transitions. However, ad ditional or larger err ors are 
likely in some cases. It was shown in lGuitou et alj Q201lh that in 
the case of transitions between the lowest-lying states, inclusion 
of other molecular states and symmetries did not signi ficantly af- 
fect th e cross-sections; this is further discussed in Be lvaev et alj 
(120121) . However, for higher-lying states, the cross-sections may 
be affected (most often increased) and this has been investi- 
gated via model estimates and included in the fluctuation fac- 
tors. In addition, there are reasons to believe that in some par- 
ticular cases the quantum chemical data are more uncertain (e.g. 
for the highest-lying states) or that the cross-sections are more 
sensitive to uncertainties in the quantum chemical data (e.g. the 
3 s 2 l S — > 3i3p 3 P° transition due to large splitting between adia- 
batic potentials at the avoided crossing). In these cases, estimates 
of the additional uncertainty are made and included in the fluc- 
tuation factors. 

The fluctuation factors would be expected to have some vari- 
ation with temperature since expected errors in the cross sections 
vary with collision energy. Generally, we expect larger errors at 
lower temperatures where the near-threshold cross-sections have 
the largest contributions. However, we have provided a single 



maximum value for all temperatures, since our estimated fluctu- 
ation factors are generally small and their variation with temper- 
ature would be even smaller, and it is not clear that any estimates 
at this level would be meaningful. 

It is worth commenting that the fluctuation factors here 
are significan t ly lo wer than those estimated for Na+H in 
Bark lem et al.l (2010). The larger values for Na+H were predom- 
inantly due to larger disagreement found between results from 
different input quantum chemical data and results from differ- 
ent quantum scattering codes. We have found no evidence for 
such large uncertainties in Mg+H, noting that the data and codes 
used here are all quite recent, while those compared in Na+H 
included older data and codes and a wider range of methods. 



3. Discussion 

The rate coefficients for Mg+H collisions are plotted against 
transition energy, AE, in Fig. Q] Thes e plots are similar t o 
those presented for Li+H and Na+H in iBarklem et alj d2011l) . 
and include rate coefficients from quantum scattering calcula- 
tions (upper panel) and those predicted by the commonly used 
Drawin formula (lower panel). Rate coefficients for charge trans- 
fer processes (ion-pair production) are shown in addition to 
those for excitation. A fit to the data for excitation processes 
is shown (dotted line ), and the fits for Li+H and Na+H from 
Bar klem et al.l (l201ll) are shown for comparison (dashed lines). 
We emphasise that the actual data have a large scatter around 
these fits, and their purpose is to indicate the general behaviour, 
not to provide fits for use in models. 

The most notable aspect of the Drawin formula results is that 
of the 21 excitation transitions considered, only 5 are optically 
allowed and can be calculated with the Drawin formula. This is 
due to the fact that Mg has two valence electrons leading to the 
presence of singlet and triplet spin terms, between which opti- 
cal transitions are forbidden. As seen in Li and Na, there is a 
tendency for the Drawin results to be generally larger than those 
from quantum scattering calculations by a few orders of magni- 
tude with a significant variation, here ranging from zero to four 
orders of magnitude larger. 

Some general features of the rate coefficients for excitation 
from quantum mechanical calculations are very similar to those 
seen in Li and Na. As expected, there is a correlation with the 
transition energy, although with a scatter of a few orders of mag- 
nitude, since the threshold energy is a main parameter in deter- 
mining the rate coefficient at low temperatures. There are also 
strong secondary correlations among transitions from a specific 
initial state. 

On the other hand, comparison of the general magnitude of 
the rate coefficients finds significant differences compared to Li 
and Na. We see that at the same AE, especially for transitions 
with large AE, the rate coefficients for Mg are generally larger 
than those for Li and Na, though converging at small AE, The 
explanation for this involves a number of factors, but the main 
reasons seem to be three-fold, and can be illustrated by consid- 
ering the transition from the ground to the first excited state. 

First, due to the different atomic structure of Mg, there is a 
tendency for analogous transitions to have larger AE in Mg than 
in Li and Na; the ground and first excited states are separated by 
2.7 eV in Mg rather than around 2 eV in Li and Na. Second, the 
larger ionisation potential of Mg compared to Li and Na means 
that avoided ionic crossings generally occur at shorter intemu- 
clear distances and therefore the radial coupling between states 
is stronger. Though the crossing occuring at smaller intern uclear 
distance generally decreases the cross-section, the increase in 
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Table 1. Rate coefficients (crv), in units of cm s , for selected temperatures in the range T - 500-8000 K, for the exci- 
tation and deexcitation processes Mg(3s nl 2s+1 L)+H(ls)^>Mg(3s n'l' 2S+1 L')+H(ls), or where indicated ion-pair production 
Mg(3s nl 2s+1 L)+H(ls)^Mg + (3i 2 S)+H~ and mutual neutralisation Mg + (3s 2 S)+R- -^ Mg(3s n'l' 2S ' +1 L')+H(ls). 



Initial 














final state 














state 


3* 21 


S 


3s3p 3 


po 


3.v3p ' 


po 


3s4s 3 S 


3i4.v 


'S 


3s3d 


'D 


3.v4p 3 P° 


Mg++IT 
















500K 
















3s 1 l S 


— 




1.18E- 


-41 


1.92E- 


-58 


5.29E-65 


548E- 


-68 


3.76E- 


-71 


4.76E-73 


8.68E- 


-83 


3s3p 3 P C 


2.95E- 


-15 


— 




7.84E- 


-28 


8.81E-35 


5.56E- 


-38 


2.99E- 


-41 


2.81E-43 


1.37E- 


-53 


3s3p 'P° 


1.10E- 


-14 


1.79E- 


-10 


— 




1.83E-17 


9.51E- 


-21 


1.38E- 


-24 


2.05E-26 


3.07E- 


-35 


3s4s 3 S 


9.71E- 


-14 


6.45E- 


-10 


5.85E- 


-10 


— 


4.33E- 


-13 


4.98E- 


-17 


2.52E-18 


6.05E- 


-27 


3i4.v'S 


2.50E- 


-13 


1.01E- 


-09 


7.57E- 


-10 


1.08E-09 


— 




1.31E- 


-12 


1.64E-15 


1.99E- 


-22 


3s3d 'D 


7.20E- 


-14 


2.28E- 


-10 


4.59E- 


-11 


5.19E-11 


549E- 


-10 


— 




4.27E-12 


341E- 


-20 


3s4p 3 P C 


3.19E- 


-14 


7.52E- 


-11 


241E- 


-11 


9.22E-11 


241E- 


-11 


1.50E- 


-10 


— 


3.29E- 


-19 


H"+Mg + 


245E- 


-13 


1.55E- 


-10 


1.51E- 


-09 


9.31E-09 
2000K 


1.23E- 


-07 


5.03E- 


-08 


1.38E-08 


— 




3/ »S 


— 




3.49E- 


-21 


248E- 


-25 


2.20E-26 


4.07E- 


-27 


4.72E- 


-28 


1.96E-28 


9.23E- 


-31 


3s3p 3 P° 


2.67E- 


-15 


— 




2.57E- 


-15 


1.04E-16 


1.20E- 


-17 


1.05E- 


-18 


344E-19 


6.15E- 


-22 


3s3p 'P° 


9.46E- 


-15 


1.28E- 


-10 


— 




6.58E-12 


4.97E- 


-13 


2.09E- 


-14 


8.84E-15 


1.30E- 


-16 


3s4.s 3 S 


6.33E- 


-14 


3.89E- 


-10 


4.95E- 


-10 


— 


6.26E- 


-11 


1.89E- 


-12 


1.98E-12 


6.68E- 


-14 


3i4.v'S 


1.88E- 


-13 


7.23E- 


-10 


6.02E- 


-10 


1.01E-09 


— 




6.27E- 


-10 


2.69E-11 


1.52E- 


-11 


3s3d 'D 


2.95E- 


-14 


8.61E- 


-11 


342E- 


-11 


4.11E-11 


848E- 


-10 


— 




2.03E-10 


8.16E- 


-12 


3s4p 3 P C 


1.92E- 


-14 


4.41E- 


-11 


2.27E- 


-11 


6.76E-11 


5.70E- 


-11 


3.18E- 


-10 


— 


3.61E- 


-12 


H"+Mg + 


2.13E- 


-13 


1.85E- 


-10 


7.88E- 


-10 


5.36E-09 
4000K 


7.60E- 


-08 


3.01E- 


-08 


8.50E-09 


— 




3s 2 l S 


— 




1.67E- 


-17 


9.32E- 


-20 


5.37E-20 


2.14E- 


-20 


6.31E- 


-21 


4.69E-21 


5.05E- 


-22 


3s3p 3 P° 


4.87E- 


-15 


— 




2.76E- 


-13 


7.95E-14 


2.07E- 


-14 


4.35E- 


-15 


2.68E-15 


147E- 


-16 


3s3p 'P 


1.05E- 


-14 


1.07E- 


-10 


— 




5.21E-11 


7.88E- 


-12 


9.96E- 


-13 


6.66E-13 


1.84E- 


-13 


3s4s 3 S 


5.26E- 


-14 


2.67E- 


-10 


4.52E- 


-10 


— 


1.38E- 


-10 


1.18E- 


-11 


1.87E-11 


9.14E- 


-12 


3i4.v'S 


1.46E- 


-13 


4.83E- 


-10 


4.75E- 


-10 


9.56E-10 


— 




1.42E- 


-09 


1.52E-10 


8.64E- 


-10 


3s3d 'D 


2.23E- 


-14 


5.28E- 


-11 


3.12E- 


-11 


4.28E-11 


741E- 


-10 


— 




4.62E-10 


1.73E- 


-10 


3s4p 3 P C 


1.55E- 


-14 


3.03E- 


-11 


1.95E- 


-11 


6.3 IE- 11 


7.38E- 


-11 


4.31E- 


-10 


— 


4.59E- 


-11 


H"+Mg + 


242E- 


-13 


242E- 


-10 


7.84E- 


-10 


448E-09 
6000K 


6.10E- 


-08 


2.35E- 


-08 


6.67E-09 


— 




3/ l S 


— 




4.52E- 


-16 


8.22E- 


-18 


7.89E-18 


3.96E- 


-18 


1.64E- 


-18 


1.34E-18 


4.56E- 


-19 


3s3p 3 P° 


9.56E- 


-15 


— 




1.29E- 


-12 


6.74E-13 


2.21E- 


-13 


6.57E- 


-14 


4.64E-14 


9.22E- 


-15 


3s3p 'P 


1.33E- 


-14 


9.86E- 


-11 


— 




1.10E-10 


1.89E- 


-11 


3.82E- 


-12 


2.70E-12 


2.34E- 


-12 


3s4s 3 S 


5.39E- 


-14 


2.18E- 


-10 


4.63E- 


-10 


— 


1.81E- 


-10 


2.66E- 


-11 


4.03E-11 


4.87E- 


-11 


3i4.v'S 


1.42E- 


-13 


3.76E- 


-10 


4.19E- 


-10 


9.50E-10 


— 




1.73E- 


-09 


2.83E-10 


3.20E- 


-09 


3s3d 'D 


2.23E- 


-14 


4.22E- 


-11 


3.20E- 


-11 


5.29E-11 


6.56E- 


-10 


— 




6.35E-10 


4.60E- 


-10 


3s4p 3 P C 


1.43E- 


-14 


2.34E- 


-11 


1.78E- 


-11 


6.29E-11 


840E- 


-11 


4.98E- 


-10 


— 


1.02E- 


-10 


H"+Mg + 


2.80E- 


-13 


2.67E- 


-10 


8.88E- 


-10 


4.37E-09 
8000K 


547E- 


-08 


2.08E- 


-08 


5.89E-09 


— 




3/ l S 


— 




3.14E- 


-15 


9.03E- 


-17 


1.02E-16 


6.16E- 


-17 


2.80E- 


-17 


2.33E-17 


145E- 


-17 


3s3p 3 P° 


1.79E- 


-14 


— 




2.79E- 


-12 


1.91E-12 


7.16E- 


-13 


249E- 


-13 


1.79E-13 


7.26E- 


-14 


3s3p 'P 


1.75E- 


-14 


9.51E- 


-11 


— 




1.71E-10 


2.97E- 


-11 


8.02E- 


-12 


5.53E-12 


8.66E- 


-12 


3s4s 3 S 


5.81E- 


-14 


1.92E- 


-10 


5.02E- 


-10 


— 


2.14E- 


-10 


4.41E- 


-11 


6.04E-11 


1.17E- 


-10 


3s4.v'S 


1.61E- 


-13 


3.28E- 


-10 


3.99E- 


-10 


9.77E-10 


— 




1.86E- 


-09 


3.93E-10 


6.08E- 


-09 


3s3d 'D 


2.36E- 


-14 


3.68E- 


-11 


348E- 


-11 


6.50E-11 


6.00E- 


-10 


— 




7.61E-10 


7.40E- 


-10 


3s4p 3 P C 


1.41E- 


-14 


1.90E- 


-11 


1.73E- 


-11 


641E-11 


9.13E- 


-11 


548E- 


-10 


— 


1.51E- 


-10 


H-+Mg+ 


3.18E- 


-13 


2.79E- 


-10 


9.79E- 


-10 


448E-09 


5. HE- 


-08 


1.93E- 


-08 


547E-09 


— 





Table 2. Estimated values of the "fluctuation factors" (a measure of the uncertainty) in the rate coefficients, given as factors of the 
rate coefficients in Table Q] Numbers above the main diagonal are the maximum values of the fluctuation factor for each transition, 
while numbers below the main diagonal are the minimum values for the same transition. 



Initial/final 








initial/final state 








state 


3s 2 'S 


3s3p 3 P° 


3s3p 'P° 


3s4s 3 S 


3s4.v'S 


3s3d 'D 


3s4p 3 P° 


Mg + +1T 


3i 21 S 


— 


1.5 


1.2 


1.2 


1.2 


1.2 


1.3 


1.2 


3s3p 3 P° 


0.6 


— 


1.2 


1.2 


1.2 


1.2 


1.3 


1.2 


3s3p 'P° 


0.8 


0.8 


— 


1.4 


1.2 


2 


1.3 


2 


3s4s 3 S 


0.8 


0.8 


0.8 


— 


1.2 


2 


1.3 


1.2 


3s4s l S 


0.8 


0.8 
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Fig. 1. Rate coefficients at 6000 K for excitation and ion-pair 
production processes in Mg+H collisions, plotted against the 
transition energy AE. Different symbols are used to denote the 
initial state of the transition following the key given in the up- 
per panel. The upper panel shows the results fr om this work 
based on detailed quantum scattering calculations (Be lvaevet alJ 
2012). Symbols inside circles refer to ion-pair production pro- 
cesses (IPP); i.e. the final state is the ionic state Mg + (3s 2 S)+H~. 
The dotted line shows a linear fit to the data for excitation pro- 
cesses (i.e. excluding ion-pair production), which is repeated in 
the lower panel to aid comparison. The dashed lines show the 
fits to the corresponding data for Li and Na from Barklem et al. 
(2011). The lower panel shows the results of the Drawin for- 
mula, and the inset shows the ratio with the quantum scattering 
calculations. Optically forbidden transitions, where the Drawin 
formula is not applicable, are shown in the bottom of the panel. 



coupling more than compensates for this. In the case of the ra- 
dial coupling between the ground and first excited state of MgH, 
the coupling peaks at around an inter nuclear distance o f 5 a.u. 
with a magnitude of around 0.4 a.u. dGuitou et al.ll2010» . while 
Li and Na both have peaks at around 7 a.u. with magnitudes of 
around 0.25 a.u. (Gadea & Boutalib 1993; Belvaev & Barklem 
2003: IM6 et aflll985t iBelvaev et al.ll2Q10h . This latter effect in 
particular leads to generally larger cross-sections between low- 
lying states. Finally, the statistical weights of the entrance chan- 
nels in the important molecular symmetry ( 2 £ + in Mg, '2 + in 
Li and Na) are generally larger in Mg than in Li and Na. For 
the case of the initial channel being the ground state, in Mg the 
2 S + symmetry is the only state available, and thus the statistical 



weight is unity. However, in Li and Na, the statistical weight of 
the initial channel is only 0.25. These factors combined lead to a 
cross-section for excitation from the ground to first excited state 
near the threshold being some two orders of magnitude larger 
for Mg than for Li and Na. This effect is reduced in the rate co- 
efficient due to the larger AE, but still leads to a rate coefficient 
roughly an order of magnitude larger. 

One can also see in Fig. Q] that the rates for charge transfer 
processes are relatively large in some cases compared to those 
for excitation, especially from the As 'S state. As in Li and Na, 
the first excited S -state again provides the largest charge trans- 
fer cross-section and rate coefficient. The rate coefficient for 
this process at 6000 K is 3.20 x 10~ 9 cm 3 s and is slightly 
larger, though of the same order of magnitude, than those from 
the first excited S -states in Li and Na, 1.72 x 10~ 9 cm 3 s -1 and 
2.03 x 10~ 9 cm 3 s -1 , respectively. The reasons that the first ex- 
cited S -state generally provides the largest rate coefficients for 
ion-pair production are that, in these three cases, the first ex- 
cited states lead to ionic crossings at intermediate internuclear 
distances where the transition probability becomes optimal (for 
detail s see discu ssion under §111. C. 2, "The second mechanism", 
in IBelvaev et alj|2012l) and that S -states lead to the largest sta- 
tistical weights for the initial channels. 

As mentioned, in non-LTE calculations for Li and Na, ex- 
citation processes were found to be less important than charge 
transfer processes; the only process with any appreciable effect 
was the charge transfer involving the first excited S -state. It is 
notoriously difficult to predict effects in non-LTE models, and 
this case is complicated by the fact that Mg, unlike Li and Na, 
has different spin terms, singlet and triplet, with large collisional 
rates which are only weakly radiatively coupled. This fact, to- 
gether with the generally higher rates for excitation processes, 
might lead one to expect that collisional excitation could be im- 
portant. However, any conclusion must await detailed non-LTE 
calculations. 
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